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ABSTRACT 

Some tails of ram-pressure stripped galaxies are detected in H I, some in Ha, and some in X-ray 
(but never all three so far). We use numerical simulations to probe the conditions for the production 
of X-ray bright tails, demonstrating that the primary requirement is a high pressure intracluster 
medium (ICM). This is because the stripped tail is mostly in pressure equihbrium with the ICM, but 
mixing leaves it with densities and temperatures intermediate between the cold gas in the disk and 
the hot ICM. Given a high enough ICM pressure, this mixed gas lies in the X-ray bright region of the 
phase diagram. We compare the simulations to observations of the ram pressure stripped tail of ESO 
137-001, showing excellent agreement in the total measured X-ray and Ha emission and non- flaring 
morphology of the tail, and consistent H I measurements. Using these comparisons we constrain the 
level of mixing and efficiency of heat conduction in the intracluster medium (ICM). 

Subject headings: galaxies: clusters, galaxies: interactions, methods: N-body simulations 



1. INTRODUCTION 

Ram pressure (and related processes) by the intraclus- 
ter medium (ICM) can remove a galaxy's gas (Gunn 
& Gott 1972). This process has been observed in var- 
ious stages; for example, Vollmer (2009) separates Virgo 
galaxies into pre-peak, peak, and post-peak ram pres- 
sure groups. The amount of time that a galaxy has been 
stripped can be estimated using the length of the observ- 
able tail and the velocity of the galaxy (e.g. Oosterloo 
& van Gorkom 2005; Sun et al. 2006). This calcula- 
tion is uncertain due to difficulties in determining the 
three dimensional galaxy velocity. Another assumption 
implicit in this calculation is that the observed tail pro- 
vides the true length of the stripped gas. In fact, tails 
have been observed in H I, Ha, and X-ray emission, al- 
though never all from the same tail (e.g. Oosterloo & 
van Gorkom 2005; Koopmann et al. 2008; Kenney et al. 
2008; Yoshida et al. 2002; Yoshida et al. 2004a,b; Sun 
et al. 2007; Sun et al. 2006; Sun et al. 2010; Machacek 
et al. 2006; Sun & Vikhlinin 2005). The lengths of tails 
observed in different wavelengths can be quite different; 
for example the H I tail of NGC 4388 is nearly three 
times as long as the observed Ha tail (Oosterloo & van 
Gorkom; Yoshida et al. 2002). Another method used to 
calculate the age of a tail is to use the estimated survival 
time of Ha, as in Gavazzi et al. (2001). However, it is 
still not clear what dictates cloud survival or even what 
conditions are necessary to produce the various types of 
emission (Ha, X-ray, and H I). Can all three types of 
emission coexist? What physical processes dominate the 
heating and mixing of stripped gas into the ICM? These 
processes include: turbulent mixing, which can generate 
intermediate temperature and density gas at constant 
pressure; shock heating, which heats the ISM; radiative 
cooling, which can lead to recompression of heated gas, 
and heat conduction, which can evaporate small clouds. 

In this work we focus on answering these questions by 
simulating gas stripping and comparing our simulated 
tail to a single observed stripped galaxy, ESO 137-001, 
which has been studied observationally in some detail. 
ESO 137-001 is in A3627, which is the closest massive 



cluster (z=0.0163, (Jradiai = 925 km s~-^ and kT = 6 
keV), similar to Coma and Perseus in mass and galaxy 
content (Sun et al. 2009 and references therein). ESO 
137-001 is a small (0.2L*; Sun et al. 2006), blue emission- 
line galaxy (Woudt et al. 2004), that is ^200 kpc from 
the center of the cluster in projection. Because its radial 
velocity is close to the average velocity of A3627 (Woudt 
et al. 2004; Woudt et al. 2008), most of its motion is 
likely in the plane of the sky, and therefore the stripping 
process is seen edge-on. Sun et al. (2006) found a ^^70 
kpc X-ray tail pointing away from the cluster center us- 
ing Chandra and XMM-Newton data. Sun et al. (2007) 
then discovered a 40 kpc Ha tail with over 30 emission- 
line regions extending through the length of the Ha tail, 
and concluded that the emission-line regions are giant 
H II regions. In a recent follow-up paper. Sun et al. 
(2009) used deep Chandra data and Gemini spectra to 
characterize the X-ray tail and H II regions in detail. 
They found a narrower secondary X-ray tail with a sim- 
ilar length. They also confirmed that 33 emission- line 
regions are H II regions, with the furthest seven regions 
beyond the tidally-truncated halo of 15 kpc that is cal- 
culated in Sun et al. (2007) using simulations by Gnedin 
(2003). In addition to these distinct H II regions, they 
find diffuse Ha emission. 

Vollmer et al. (2001) searched for H I in A3627, and 
did not detect any H I in or around ESO 137-001 with 
a limiting column density of 2 x 10^^ cm~^ and a reso- 
lution of 15". In fact, of the ~80 galaxies identified by 
Woudt et al. (1998) in their search region, Vollmer et 
al. (2001) detected only 2 in H I, finding that the H I 
detection rate in A3627 is similar to that in Coma. 

Sivanandam et al. (2009) observed ESO 137-001 with 
IRAC and IRS on Spitzer. The IRS data extended to 
20 kpc from the galaxy along the X-ray tail, and warm 
(~160 K) molecular Hydrogen was detected throughout 
the length of the observed region. The observed region 
contains ~2.5 x 10^ Mq warm H2 gas. They also iden- 
tify star-forming regions using 8 jam data, which coincide 
with Ha emitting regions. 

There has been a substantial amount of theoretical 
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TABLE 1 

Galaxy Stellar and Dark Matter Constants 



Variable 


Value 




1 X IQii M© 




4 kpc 


b* 


0.25 kpc 


^bulge 


1 X IQio Mq 


^bulge 


0.4 kpc 


^DM 


23 kpc 


PDM 


3.8 X 10-25 g 



work investigating ram pressure stripping in general (e.g. 
Schulz & Struck 2001; Quilis, Bower & Moore 2000; 
Roediger & Briiggen 2008, Kronberger et al. 2008; 
Kapferer et al. 2009) - see Tonnesen & Bryan (2009, 
2010; hereafter TB09 and TBIO) for a more detailed 
discussion. There have also been simulations designed 
to predict or interpret observational characteristics of 
ram pressure stripped tails and the remaining disks (e.g. 
Vollmer et al. 2005, 2006, 2008), but detailed, quanti- 
tative predictions of all three observational probes have 
been missing to date (H I, diffuse Hce, and X-ray emis- 
sion) . 

In our previous work (TBIO), we ran a set of high reso- 
lution simulations (about 38 pc resolution, which is small 
enough to marginally resolve giant molecular clouds) to 
understand how a multiphase ISM could affect the sur- 
vival and structure of ram pressure stripped gas. We 
focused on how density fluctuations that are observed 
in the multiphase ISM of galaxies can affect gas tails. 
Including radiative cooling allowed us to estimate the 
density of and emission from H I, Ha, and X-ray gas 
separately. We found that both the morphology and ve- 
locity structure of our tails agreed with observations of 
long gas tails (e.g. Oosterloo & van Gorkom 2005). Our 
simulations also resulted in observable amounts of H I 
and Ha emission. However, the X-ray tail had a low sur- 
face brightness, which we attributed to the low pressure 
of the surrounding ICM. In this paper we will use the 
same method as in TBIO, but have chosen ICM param- 
eters comparable to the ICM around ESQ 137-001. By 
focusing on the level of agreement between our simula- 
tions and the observations of ESO 137-001 we will be able 
to discuss the importance of physical mechanisms such 
as heat conduction, as well as predict the conditions un- 
der which H I, Ha, and X-ray emission are produced in 
stripped tails. 

The paper is structured as follows. After a brief in- 
troduction to our methodology, we provide the charac- 
teristics of our simulations and our method of producing 
simulated observations (§2). We then (§3) present our re- 
sults, specifically focusing on the comparison with ESO 
137-001. In §4 we discuss the broader implications of 
our simulation, and discuss our choice of radiative cool- 
ing floor and resolution in §5.1-2. Finally, we conclude 
in §6 with a summary of our results and predictions for 
observers. 

2. METHODOLOGY 

We use the adaptive mesh refinement (AMR) code 
Enzo (Bryan 1999; Norman & Bryan 1999; O'Shea et 
al. 2004). Our simulated region is 311 kpc on a side 



TABLE 2 
Gas Disk Constants 



Variable 


Value 


MgaS 


1 X IQio M0 


^gas 


7 kpc 




0.4 kpc 



with a root grid resolution of 128^ cells. We allow an 
additional 6 levels of refinement, for a smallest cell size 
of 38 pc. We refine our simulation based on the local gas 
mass, such that a cell was flagged for refinement when- 
ever it contained more than about 2 x 10^ Mq. We found 
that this refined most of the galactic disk to 38 pc resolu- 
tion by the time the wind hit; dense clumps in the wake 
were also refined to 38 pc resolution, while more diffuse 
components had lower resolution. 

The simulation includes radiative cooling using the 
Sarazin & White (1987) cooling curve extended to low 
temperatures as described in Tasker & Bryan (2006). To 
mimic effects that we do not model directly (such as tur- 
bulence on scales below the grid scale, UV heating, mag- 
netic field support, or cosmic rays), we cut off the cool- 
ing curve at a minimum temperature Tmin so that the 
cooling rate is zero below this temperature. In the sim- 
ulations described here we use either Tmin = 8000 K, or 
Tmin = 300 K. Both allow gas to cool below the threshold 
for neutral Hydrogen formation. In TB09, we found that 
the minimum temperature affected the range of masses 
and sizes of clouds forming in the disk, and the distribu- 
tion of clouds throughout the disk, although the range 
of gas densities was very similar. The Tmin = 300 K case 
resulted in a more fragmented disk whose small clouds 
took longer to strip. Therefore the timescales for gas 
stripping differed, although the total amount of gas lost 
was similar. In the simulations presented in this paper, 
we find that the cooling floor affects the amount of frag- 
mentation in the disks, but not the mass-distribution of 
gas densities. The larger radii of the remaining gas disks 
in the two Tmin = 300 K runs are due to the survival of 
dense clouds in the outer disk (Figures 1, 4, and 5). 

In TBIO, we found that the structure of the wake also 
depends somewhat on Tmin- Since the publication of 
TBIO, we found and corrected an error in our implemen- 
tation of the cooling rate that resulted in a slight shift 
in the peak of the cooling curve around 10^ K, but did 
not significantly affect cooling at low and high temper- 
atures (this affected only the Tmin = 300 K simulation). 
Tests showed that this had only a small impact on the 
dynamics of the flow, and on the predicted H I and X-ray 
measures, but did strongly affect the predicted Ha emis- 
sion, which is extremely sensitive to the gas temperature 
around 10^ K. This explains the enhanced Ha emission 
in this paper for the Tmin = 300 K run as compared to 
TBIO. 

Our galaxy model is the same as in TBIO and TB09, 
which used the spiral galaxy model described in Roedi- 
ger & Briiggen (2006). We list the model parameters 
for our galaxy in Tables 1 and 2. The stellar and dark 
matter components of the galaxy are static potentials. 
As in TB09 and TBIO, our galaxy position and box size 
allows us to follow gas 200 kpc above (in the wind, or 
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TABLE 3 
Runs summary 



Run 


vicM (km/s) 


PiCM (dyne/cm^) 


TicM (K) 


Pram (dyne/cm^) 


tproj (Myr) 


T80vh 


1900 


4.2 X 10-11 


8.3 X 10'^ 


11.6 X 10-11 


85 


T3vh 


1900 


4.2 X 10-11 


8.3 X 10'^ 


11.6 X 10-11 


75 


TSvl 


1413 


2.66 X 10-11 


7.3 X 10'^ 


5.29 X 10-11 


110 


Sun et al. 2010 




1.8 X 10-11 


-7 X 10^ 







z, direction) the galaxy. To identify gas that has been 
stripped from the galaxy, we also follow a passive tracer 
which is initially set to 1.0 inside the galaxy (defined as 
gas that is above the ICM density) and 10~^ outside. 
In the following analysis, we will use a minimum tracer 
fraction of 0.25 to find gas stripped from the galaxy (as 
in TBIO). 

2.1. Introduction to the Three Runs 

The galaxy initially evolves in a static, high-pressure 
ICM with p = 8.7 X 10-^^ g cm"^ and T = 9.069 x 
10^ K {p =7.57 X 10-2^ g cm-3 and T = 1.04294 x 10^ 
K for the slower wind case), to allow cool, dense gas to 
form in the galaxy (each of our three runs has about 3 x 
10^ M0 of gas with densities at or above 10~^^ g cm~^ 
when the wind hits the disk) . This naturally generates a 
multiphase ISM (see Tasker & Bryan 2006 and TB09 for 
more discussion of the ISM properties). 

After 155 Myrs, we reset the boundary conditions to 
generate a constant ICM inflow along the inner z-axis, 
which is always face-on to the galaxy. In Table 3 we show 
the details of each of the three runs. The table includes 
the ICM parameters and the time after the wind has hit 
the galaxy at which the X-ray tail is 80 kpc long (this is 
how we choose the outputs to compare to the observa- 
tions of ESQ 137-001). 'T80' indicates cooling to 8,000 
K and 'T3' indicates cooling to 300 K, while 'vh' and 'vl' 
indicate high and low velocity wind, respectively. See 
TB09 for other details regarding the general numerical 
setup, and TBIO for a discussion of the general impact of 
the cooling floor on the tail structure. In order to com- 
pare with both observations and our previous work, we 
use a face-on wind direction. 

2.2. Projections 

Enzo outputs the density and temperature of the gas 
in each cell. To transform these values into H I column 
density and Ha intensity, we used Cloudy, version 08.00 
of the code last described by Ferland et al. (1998). Us- 
ing a table of temperatures and densities, we calculated 
the hydrogen neutral fraction and Ha emissivity. This 
table is then used to calculate the observational quanti- 
ties for each cell, which are then summed to generate an 
image. This is described in detail in TBIO, but briefly, 
we included CMB radiation, the cosmic ray background, 
bremsstrahlung radiation from the ICM and the 2005 
version of the Haardt & Madau (2001) z = metagalac- 
tic continuum, as implemented by Cloudy. 

We chose to calculate the neutral fraction and Ha emis- 
sivity for a plane-parallel gas cloud of width 100 pc. We 
selected this width because it loosely corresponds to the 
cell size of most of the gas in the tails, and accounts ap- 
proximately for attenuation of the ionizing background 



radiation. If we assumed the radiative thin limit (by us- 
ing a very small cloud size in Cloudy) , it would somewhat 
decrease the amount of H I we predict, and increase the 
Ha emission for dense, low-temperature gas. We discuss 
the use of diff'erent cloud sizes in detail in TBIO, how- 
ever, in the simulations in this paper we flnd that using 
a 10 pc cloud does not signiflcantly change our Ha flux 
nor our H I column densities because of the higher densi- 
ties in our tail gas. This indicates that radiative transfer 
effects are not that important for this work. 

To create X-ray surface brightness projections, we use 
a spectral lookup table that depends on temperature and 
density, assuming a constant metallicity of 0.3 solar, as 
computed using a Raymond-Smith code (Raymond & 
Smith 1977), as updated in XSPEC (Arnaud 1996). The 
X-ray band we use is 0.5 keV to 2.0 keV, following Sun 
et al. (2006). 

3. COMPARISON TO OBSERVATIONS 

In this section, we compare our simulated stripped tail 
to observations of ESO 137-001. We choose an output 
time at which the X-ray tail length is about 80 kpc in 
order to match the observations of Sun et al. (2010); this 
time is shown in Table 3. We are using this comparison 
with the observations of ESO 137-001 to better under- 
stand the physics at work in producing X-ray, Ha, and 
H I tails, so we have not tuned our simulation specifically 
to find an exact match to the galaxy and tail of ESO 137- 
001. While our ICM parameters are similar to those of 
the ICM near ESO 137-001, they are not the same. The 
final row in Table 3 displays the ICM parameters from 
the observations of Sun et al. (2010). We have not at- 
tempted to model the exact angle between the galaxy's 
disk and orbital motion (Woudt et al. (2008) find a posi- 
tion angle of 125°), and instead use a face-on wind, and 
we are using a larger galaxy (Initially our galaxy has a 
radius of 26 kpc, and in the comparison projections the 
radius is about 15 kpc, while the 2MASS isophotal ra- 
dius for ESO 137-001 is 6.1 kpc (Skrutskie et al. 2006)). 
In the following comparisons we will take note of the dif- 
ferences between our simulations and the observations of 
Sun et al. (2006; 2007; 2010). 

We have chosen the output time at which each simula- 
tion has an X-ray tail whose length is about 80 kpc. This 
means that the length of time that each galaxy has been 
stripped in the different simulations is not the same. In 
Table 4, we list the amount of gas in the tail in three dif- 
ferent temperature ranges. These masses include all of 
the gas above 10 kpc from the disk with tracer fractions 
of at least 25%, so some of the gas included is at low 
density, as shown in Figure 2 (too low to be observable). 
Nevertheless, it can be very loosely stated that longer 
stripping times may result in more mixing of hot gas 
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Fig. 1. — X-ray surface brightness maps of our three simulated galaxies, seen side-on. From left to right, we show T3vl, T80vh, T3vh. We 
have chosen the outputs with tail lengths of 80 kpc for our comparison to the observations of Sun et al. (2006; 2007; 2010). See discussion 
in Section 3.1.1 



TABLE 4 

Amount of Stripped gas at Different Temperatures 



Run 


T < 10^ K 
10^ M0 


10^ < T < 10^ K 
10^ M0 


7 X 10^ < T < 4 xlO'^ K 
10^ Mq 


TSOvh 


1.4 


1.6 


6.0 


T3vh 


1.1 


1.2 


9.1 


T3vl 


4.2 


2.9 


5.9 



into the ICM, and in more dense gas being stripped and 
condensing into clouds in the tail. This can be seen by 
comparing TSOvh or T3vl with T3vh (of course there are 
also other differences in the simulations that influence 
the tail gas, such as cooling floor and ICM pressure). 
However, we stress that in this paper we mainly focus on 
what causes a tail to be bright in X-ray and Ha emission, 
so focus more on ICM pressure than the amount of time 
a galaxy has been stripped (although we do consider this 
timescale in our discussion of heat conduction). 

3.1. X-ray 
3.1.1. Comparison to ESQ 137-001 

We first compare the X-ray characteristics of our sim- 
ulations to observations. The X-ray surface brightness 
projections of our runs are shown in Figure 1. As in 
Sun et al. (2006; 2010), we measure the X-ray emis- 
sion between 0.5-2.0 keV. The surface brightness profiles 
along the tails are in rough agreement with that of ESO 
137-001, which is bright near the disk and has a second 
bright region about 40 kpc from the disk before becom- 
ing less luminous to the end of the tail. Our model does 
not reproduce the exact X-ray surface brightness profiles: 
we do not have a model that includes both the bright 
emission near the disk and a surface brightness decrease 
across the entire tail before the second brightness peak. 

To make the X-ray projections, we adopted a minimum 
observable surface brightness of 7.1 x 10~^ erg cm~^ 
s~^, which we estimated from the total luminosity of 



the observed tail and brightness profiles in Sun et al. 
(2010) (our cluster background is about 3.6 x 10~^ erg 
cm~^ s~^). This results in a luminosity difference of 
an order of magnitude between the lowest and highest 
surface brightness features in all of our simulated tails, 
similar to that in the observed tail of ESO 137-001. 

Our tails are narrow and show a nearly constant width 
along the entire tail, in agreement with the X-ray mor- 
phology reported in Sun et al. (2006; 2010). This is 
in contrast to simulations that do not include radiative 
cooling and produce flared tails (e.g. Roediger, Briiggen 
& Hoeft 2006) However, it is clear that we do not have 
separated X-ray tails as in the observations. This is not 
surprising given the explanation of Sun et al. (2010) 
that the two tails likely result from the stripping of two 
spiral arms. As we discuss in detail in TB09, our disks 
fragment but do not form spiral arms. 

In Table 5, we list some of the characteristics of our 
simulated X-ray tails to compare with the observed tail. 
To calculate the luminosity of the tail, we find the total 
energy emitted and subtract the background emission 
from the ICM. As the table demonstrates, the simulated 
tails are wider, more luminous, and have a higher aver- 
age (emission-weighted) temperature than the (spectro- 
scopic) temperature measured by Sun et al. (2010). We 
now address these differences. 

The tail is most likely wider because we are modeling 
a large spiral galaxy, while ESO 137-001 is thought to 
be a ~ 0.2L* galaxy, with a smaller galactic radius (Sun 
et al. 2006 and references therein). In fact, the entire 
volume of the tail in any of our three runs is nearly an 
order of magnitude larger than the tail of ESO 137-001 
(we assume that the tails are cylindrical with the heights 
and diameters denoted in Table 5). The other difference 
is that our ICM pressure is also somewhat larger than 
calculated in Sun et al. (2010) by either a factor of 1.47 
(T3vl) or 2.33 (T3vh and TSOvh). 

These differences in tail volume and ICM pressure im- 
pact the observables in two ways. First, as we discuss 
in TBIO and below, the compression of stripped gas by 
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TABLE 5 

X-RAY TAIL ATTRIBUTES 



Run or 


/ X w 


Lo.5-2keV 


a 

Lo.5-2 corrected 


T 


Observation 


(kpc X kpc) 


(10^0 erg s-1) 


(10^0 erg s-i) 


(10'^ K) 


T80vh 


80 X 26 


66.3 


4.8 


2.5 


T3vh 


80 X 30 


93.0 


5.1 


2.1 


T3vl 


80 X 30 


41.9 


3.6 


1.4 


Sun et al. 2010 


80 X 8, 80 X 7 


(8.3 d= 0.4) 


(8.3 ± 0.4) 


0.93 ± 0.05 



^The corrected luminosity multiplies the simulated X-ray lumi- 
nosity by the ratio between the simulated and observed tail volumes 
and the ratio between the simulated and observed ICM thermal 
pressures (see text). 



the ICM determines the temperature and density dis- 
tribution of gas in the tail. Therefore, a higher ICM 
pressure results in higher-density hot gas (T > 10^ K), 
which will have a higher X-ray emissivity. As we argue 
in the next section, this makes the X-ray luminosity pro- 
portional to the ICM pressure. Second, if we assume 
that the filling factor of X-ray emitting gas is the same 
in our simulations as in the tail of ESO 137-001, then 
the total luminosity is also directly proportional to the 
volume of the entire tail. The X-ray luminosity after 
applying these corrections is also shown in Table 5 (as- 
suming X-ray luminosity is directly proportional to both 
the ICM pressure and the volume of the tail). When we 
account for these differences, our measured luminosities 
are within about a factor of two of the X-ray luminosity 
of ESO 137-001. 

As shown in Table 5, there is also a factor of two dif- 
ference between the temperatures of the simulated and 
observed X-ray tails; however, this is most likely due 
to the different ways in which this quantity is deter- 
mined in simulations as compared to observations. We 
use luminosity- weighted temperatures, which tend to be 
higher than spectroscopic temperatures when there is a 
range of gas temperatures (see Mazzotta et al. 2004). 
Sun et al. (2010) discuss how the spectral fitting of the 
iron-L hump biases their temperatures low if there are 
significant emission components at A;T = 0.4 - 2 keV, 
which is certainly the case in our simulation and in the 
tail of ESO 137-001 (Sun et al. 2010; Sivanandam et al. 
2009). Therefore, the agreement is probably consider- 
ably better than indicated in Table 5. 

3.1.2. What Makes a Tail X-ray Bright? 

In this section, we examine what lights up an X-ray 
bright tail, and explain why some observed tails, like ESO 
137-001, exhibit X-ray emission, while others do not. In 
Figure 2 we show the mass-weighted distribution of den- 
sity and temperature of gas in the wake that originated 
from the galaxy from our T3vl run (left) and from the 
^min = 300 K run from TBIO using the corrected cooling 
curve (right). The plots include all of the gas located be- 
tween 10 kpc and 240 kpc above the disk that has at least 
25% of it's mass originating in the galaxy (as determined 
by the tracer fraction). We choose to highlight these two 
cases because the ICM pressure is the only difference be- 
tween the two simulations (the thermal pressures differ 
by a factor of 15). As discussed in more detail in TBIO, 
the (T > 10^ K) gas in the wake is largely in pressure 
equilibrium and so falls roughly along a line of constant 



pressure. The impact of the higher ICM pressure can be 
seen as a shift to higher density and temperature in the 
left panel. 

In red we also plot two contours of constant luminosity 
per mass (or emissivity per density). The lower contour 
is at 10~^ erg s~^ g~-^. If the minimum observable X-ray 
surface brightness is 10~^ erg s~^ cm~^, then in order for 
X-ray emission to reach this level, there must be a suffi- 
cient amount of gas at or above the red contour. More 
precisely, we need a column density of hot gas of 10^^ 
cm~^ in order to produce observable emission (this cor- 
responds to about 10^ Mq along a single line of sight 
through our 38 pc x 38 pc cells). In general, only the 
colder gas (T < 10^ K) in the tail has these high sur- 
face densities (lO^-*^ cm~^), so we also plot an emissivity 
line which only requires a surface density of 10^^ cm~^ 
in order to be observable (the upper 10~^ erg s~^ g~^ 
contour). Gas that is hot enough to emit X-rays (hotter 
than ~7 x 10^ K) has this lower surface density in the 
tail (10^0 cm-2). 

This figure shows that gas at higher densities and lower 
temperatures than the ICM (but above ^ 10^ K) will be 
emitting the most strongly in X-rays. As galactic gas 
is stripped it either cools into clouds (T < 10^ K) or is 
compressed to the ICM pressure and begins to mix with 
ICM gas. While the tail gas is at the high ICM pressure 
- but before it is completely mixed with the ICM - it 
will be X-ray bright. The mixing occurs along the line 
of constant pressure seen in Figure 2. 

The separation between cold clouds and X-ray emitting 
gas is seen clearly in the left panel of Figure 3, which 
shows a thin slice from a detailed section above the disk. 
The cold clouds are black in this figure, indicating no X- 
ray emission, while the X-ray emitting gas is not confined 
to regions close to the dense clouds. Most of the X-ray 
bright cells in our simulation have between 70% and 90% 
of their gas originating from the galaxy. In the lower left 
panel of Figure 3 we show the gas temperature, which 
when compared with the X-ray slice again highlights that 
while some mixing and radiative cooling can enhance the 
X-ray emission, but too much lowers the density of the 
hot gas to the point where the X-ray emissivity is low. 
The X-ray bright gas may have either been stripped as 
hot gas from the disk and slowly mixed, or stripped from 
the dense clouds in the tail and mixed into the ICM. 

3.2. Ha 

Next, we turn to Ha emission, shown in Figure 4. The 
minimum observable surface brightness that we adopt 
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Fig. 2. — Contour plots showing the mass in gas at different densities and temperatures for the T3vl run on the left, Tmin = 300 K (from 
TBIO) on the right. The contours are spaced by a factor of 10 in mass. For a gas cell to be included in the contour plot, it must have at 
least 25% of its mass originating from the galaxy (based on the tracer fluid). The two curves denote lines of constant X-ray luminosity per 
mass (or emissivity per density). This flgure illustrates that high ICM pressure produces an X-ray bright tail. See discussion in Section 
3.1.2 



for these maps is 2 x 10~^^ erg s~^ cm~^ arcsec"^ (see 
Sun et al. (2007) and references therein). As described 
ear her, we use Cloudy to determine the Ha emissivity 
given a gas temperature and density. 

We find, as shown in Figure 4, highly structured, long 
tails of Ha emission. Note that we do not include UV 
radiation from star formation or AGN (except from the 
metagalactic background, as described in section 2.2). 
Figure 3 shows that Ha emission peaks around the edges 
of cold clouds, which can be seen by comparing the upper 
right panel of Ha emission with the lower right panel 
showing gas surface density. 

The Tmin = 8000 K and T^in = 300 K runs show a 
similar amount of Ha emission. This is because, while 
the minimum temperature affects the temperature in the 
central regions of the clouds, it does not strongly af- 
fect the cloud edges in the simulation, whose characteris- 
tics are more determined by the interaction between the 
cloud and the ICM. We find that changing the cloud size 
parameter in our Cloudy run has only a small effect on 
the total Ha fiux (less than 10% in all three runs). This 
is because most of our emission is produced by collisional 
processes rather than photoionization, and so is mostly 
dependent on the gas temperature, not the optical depth 
to ionizing photons (recall that we do include star forma- 
tion in the simulation and so do not model HH regions 
within the tail). 

Because we do not include star formation, we cannot 
compare our Ha emission to the 33 H II regions seen by 
Sun et al. (2007). We do compare the total fiux from 
our tail to the observed diffuse Ha emission in Table 6. 
Again, the volume of our diffuse tail is much larger than 



TABLE 6 

Ua TAIL ATTRIBUTES 



Run or 


/ X w 


fHc./10-i4 


fHa/lO corrected"^ 


Observation 


(kpc X kpc) 


(erg s"-*- cm~^) 


(erg s"-*- cm~^) 


T80vh 


67 X 27 


75.8 


2.23 


T3vh 


65 X 29 


75.9 


2.00 


T3vl 


80 X 30 


97.9 


1.96 


Sun et al. (2007) 


-40 X 6 


4.4 


4.4 



^The corrected flux multiplies the simulated Ho; flux by the ratio 
between the simulated and observed tail volumes. 



that observed by Sun et al. (2007), and we show the 
corrected Ha flux by dividing by the volume ratio in the 
table. The difference in tail widths is from the difference 
in galaxy sizes, while the length of the observed tail is 
the minimum length of the diffuse emission because there 
is a bright star in the fleld (Sun et al. 2007). When we 
take the different volumes of the tail into account, our 
simulated Ha flux is less than that observed from the tail 
of ESO 137-001, but in all three cases, the simulated Ha 
flux is within a factor of 2.5 of the observed flux. 

There are a number of possible explanations for our 
slightly lower predicted Ha flux, including the possibility 
that some unresolved H II regions are counted as diffuse 
flux in the observations. Other heating sources such as 
thermal conduction — an effect we do not include in 
the simulations — could give rise to more Ha flux than 
we see in our simulations. We also flnd that numerical 
resolution plays a role in the total Ha emission, as we 
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Fig. 3. — Projections of a very thin (0.3 kpc) slice from the T80vh run (images are 16 kpc x 19 kpc). The upper left panel shows X-ray 
surface brightness (red brightest), which demonstrates that X-ray emitting gas is associated with diffuse, mixed gas, but not necessarily 
material recently stripped from dense clouds in the tail. The lower left panel (temperature) demonstrates that the brightest X-ray emitting 
gas is associated with intermediate temperatures. Similarly, the upper right panel shows that Ho; emission is produced primarily at the 
edges of dense clouds. The lower right panel (gas surface density) demonstrates that the brightest Ha is produced in high surface density 
clouds. 



will discuss in Section 5.2. 

3.3. HI 

In this section we consider H I column density, com- 
paring projections of our simulations to observations. In 
Figure 5 we show the H I column density at the best 
resolution of our simulation. As discussed in Section 2.2, 
we use Cloudy to determine the neutral fraction given a 
temperature and density (with the assumed metagalac- 
tic UV flux) and then apply that value in our projection 
routine. 

In Figure 6 we only show a projection comparable to 
the observations performed by Vollmer et al. (2001), us- 



ing a resolution of 30" and a minimum column density 
of 2 X 10^^ cm~^. We also correct for the difference 
in disk sizes by dividing our column densities by a fac- 
tor of 3 (the ratio of X-ray tail diameters). Only the 
tail of T3vl would have been observed by Vollmer et al. 
(2001). This is because in the slower velocity case there is 
a smaller velocity difference between the clouds and the 
ICM wind and the ICM is slighty less dense (83% of the 
ICM density in the higher velocity cases), which results 
in less cloud ablation and more high density gas (both 
a slightly higher maximum density of gas and slightly 
more gas at any particular density above n ~ 1 cm~^). 
This difference between the non-detection by Vollmer et 
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Fig. 4. — Ha intensity projections. These can be compared with the diffuse emission observed by Sun et al. (2007). From left to right: 
T3vl, T80vh, T3vh. See Section 3.2 for discussion. 
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Fig. 5. — H I column density projections at the maximum resolution of our simulation (38 pc). From left to right: T3vl, TSOvh, T3vh. 
See Section 3.3 for discussion. 



al. (2001c) and our prediction that T3vl would be de- 
tected implies that ESO 137-001 is most likely moving 
more quickly through the ICM than 1413 km s~^. The 
H I column densities in T3vh and T80vh are within a fac- 
tor of 3 and 1.3 respectively below the maximum column 
density as determined by the non-detection of Vollmer et 
al. (2001). 

Figure 3 clearly shows that bright Ha emission is pro- 
duced only at the edges of dense neutral clouds. This 
spatial correspondence between neutral and Ha emitting 
gas agrees with the observations of Sivanandam et al. 
(2009), who found molecular hydrogen in the tail to the 
farthest distance they could search-20 kpc. We predict 
that with a deeper observation H I will be observed to 
at least 40 kpc (the length of the observed Ha tail), and 
likely even farther (our shortest H I tail is 65 kpc in 
T3vh), unless heat conduction is quite efficient. 



4. DISCUSSION 

4.1. Which Stripped Galaxies will have X-ray Tails? 

In this section, we comment on the likelihood of finding 
more X-ray tails. In order for a tail to have observable X- 
ray emission, it must be in a high pressure ICM. We can 
estimate the pressure necessary to produce bright X-ray 
emission simply by shifting our mass-weighted distribu- 
tion of gas until the lowest mass contour lies along the 
upper 10~^ erg s~^ g~^ contour. We choose this contour 
because in every simulation we have run, projections of 
surface density have shown that most of the tail has a 
column density of at least 10^^ cm~^. Gas in the tail 
with higher surface densities tends to be in cold clouds 
(T < 10^ K) (see TBIO), while X-rays are emitted by 
hot (T > 10^ K), diffuse gas with lower column densi- 
ties (Figure 2). This method results in a minimum ICM 
pressure of 9 x 10~^^ erg cm~^. 
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TABLE 7 

Where in clusters will tails be X-ray bright? 



Fig. 6. — H I column density projection smoothed to the resolu- 
tion of the observations by Vollmer et al. (2001) using ATCA. The 
only run with observable H I in the tail using 30" resolution and a 
mmimum column density of 2 x 10^0 is T3vl. See Section 3.3 for 
discussion. 



In Table 7 we show the cluster radius at which the ICM 
pressure falls below our minimum value. In general, we 
calculated the cluster radius using a /3-model and con- 
stant temperature from the papers cited in column 3. 
There are three exceptions to this. First, the X-ray inten- 
sity contours in Figure 17 of Wang et al. (2004) denote 
the highest density ICM region (and are not spherically 
symmetric), which is the only region of A2125 with a 
high enough density to produce our minimum ICM pres- 
sure for X-ray bright tails. Second, we find the radius of 
Perseus with our minimum ICM pressure directly from 
Figure 9 in Ettori et al. (1998). Finally, we use the en- 
tropy measurements of Cavagnolo et al. (2009) to find 
that as close as 100 kpc from M87 the ICM pressure is 
below 9 X 10~^^ erg cm~^. Note that these are first 
estimates generally using spherically symmetric cluster 
profiles, which oversimplify the structure of the ICM. 
We can now comment on whether observed X-ray tails 
are in high pressure ICMs (defined here as at or above 9 
X 10"^^ dyn cm-^). The two X-ray tails in A3627, ESQ 
137-001 and ESQ 137-002, are both within the high pres- 
sure ICM (Sun et al. 2007; Sun et al. 2010). However, 
Sun et al. (2010) find that using their spectroscopic tem- 
peratures, the tail gas is over-pressured with respect to 
the ICM. This might indicate that the spectroscopic tem- 
perature does not correctly model the tail temperature. 
We find that turbulent pressure is not strong enough to 
greatly increase the gas pressure in the tail (Figure 2 and 
TBIO). 

CI 53 is near the X-ray peak of A2125 (Wang et al. 
2004). UGC 6697 in A1367 is at least 450 kpc from 
the X-ray peak of the cluster, outside of our calculated 
high-pressure radius. However, it is near an infalling sub- 
cluster and in a higher-density ICM than predicted using 
the /3-model centered on the X-ray peak of A1367. Using 
the fit to the subcluster by Donnelly et al. (1998), Sun 
& Vikhlinin (2005) calculated the surrounding pressure 
to be 7.9 X 10~^^ dyn cm~^, very close to our minimum 



Cluster 


Outer radius (kpc) 


Reference 


A3627 


250 h-,' 


Bohringer et al. (1996); 




Sun et al. (2010) 


A2125 


28.5 - 89 


Wang et al. (2004) 


A1367 


274 h-^ 


Mohr et al. (1999) 


CL 0024+0016 


250 h~Q 


Zhang et al. (2005) 


Coma 


613 h-^ 


Briel et al. (1992) 


Perseus 


2000 h~Q 


Ettori et al. (1998) 


Virgo (M87) 


< 100 


Cavagnolo et al. (2009) 



pressure for observable X-ray emission. 

4.2. Efficiency of Heat Conduction 

Heat conduction, which we do not include in these 
simulations, could be important for the survival of cool 
clouds in the ICM and Ha emission; however, it can be 
suppressed by magnetic fields (e.g. Vollmer et al. 2001). 
If it is an efficient way to transport heat from the ICM 
to cold, stripped gas, then the survival time of H I clouds 
would be less than predicted in this paper, and the length 
of the tails would be shorter. We can estimate the effi- 
ciency of heat conduction by comparing an analytic cal- 
culation of the evaporation time of the most dense clouds 
in our simulations to the length of time we expect clouds 
have survived in order to produce the observations of 
ESQ 137-001. 

First, we estimate the evaporation time for a typical 
cloud if heat conduction is not suppressed, and we will 
define this calculated evaporation time as the time for 
cloud evaporation if heat conduction is 100% efficient. 
We follow Cowie & McKee (1977), as in Vollmer et al. 
(2001). In our clouds, we find that the mean free path 
for ions is comparable to or greater than the temperature 
scale length, so we need to use the saturated heat fiux 
equations. 

Evaporation time is proportional to / 

^iCM^^ic^v/^' where / is the conduction efficiency rel- 
ative to Spitzer. Solving for the evaporation time for 
each of our three runs, we find that using a cloud radius 
of 100 pc, and the maximum gas density in our tails at 
the time of our comparison, the evaporation times of our 
three runs are - 4 Myr for T80vh, - 8 Myr for T3vh, 
and ~ 10 Myr for T3vl. We choose to use the highest 
density found in our tails because we want to calculate 
the longest plausible evaporation time in order to find 
the most conservative estimates for the maximum effi- 
ciency of heat conduction. The Tmin = 300 K runs have 
a maximum density in their tails of 2 x 10~^^ g cm~^ 
and TSOvh has a maximum density of 10~^^ g cm~^. 

We have shown that diffuse Ha emission directly traces 
neutral clouds (Figure 3), which means that the dense 
clouds must not be entirely evaporated before they reach 
40 kpc above the disk. We are also able to measure the 
velocity of gas in our tails (see TBIO for details), so can 
estimate the time it would take for a cloud to reach 40 
kpc above the disk. We use generous estimates of 1200 
km s~^ for the high velocity cases and 900 km s~^ for the 
lower velocity run in order to calculate the shortest time 
it would take these clouds to reach 40 kpc above the disk 
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(which win maximize the value we find for the efficiency 
of heat conduction). Therefore, the efficiency of heat 
conduction must be less than 24% (T3vl and T3vh), or 
12%(T80vh). 

Once there is a deep observation of the H I tail, we 
win be able to compare the length of ESO 137-001 's tail 
to the lengths of our simulated tails to ffnd a minimum 
efficiency of heat conduction. To do this we will use a 
similar argument to the one above; namely, compare the 
amount of time a cloud takes to reach its height above 
the disk to the evaporation time, assuming that heat 
conduction is responsible for destroying the clouds. It is 
important to consider our simulations in this calculation, 
because they show that even without heat conduction, 
H I may not extend along the entire length of the X-ray 
tail. 

The simulations also show the limitations of using a 
simple velocity and distance argument to calculate the 
length of time a galaxy has been stripped-using 900 km 
and the length of the H I tail we find that the galaxy 
in T3vl has been stripped for 87 Myr, while the actual 
time after the wind has hit the galaxy is 110 Myr. The 
high velocity cases have an even larger discrepancy be- 
tween the actual stripping time and the time calculated 
using the H I tail length. 

Despite these caveats, we can make a very rough first 
estimate using Figure 6 (T3vl) and the non-detection re- 
sult of Vollmer et al. (2001). The H I in our smoothed 
projection begins about 50 kpc above the disk. If the 
non-detection of Vollmer et al. (2001) means that there 
are no neutral clouds at this height above the disk, then 
we can calculate a minimum efficiency for heat conduc- 
tion, which is 18.5%. However, as we discussed above, 
we expect that a deeper observation will find H I beyond 
this height above the disk. 

5. NUMERICAL CONSIDERATIONS 

5.1. Radiative Cooling Floor 

As noted earlier, we used two different values for our 
radiative cooling ffoor (8000 K and 300 K) in order to 
explore in a simple way the potential impact of processes 
which we do not include in the simulation. We discuss the 
physical motivations for the two cooling ffoor s in TBIO. 
We ffnd that the change in the minimum temperature 
makes only a relatively minor change in the morphology 
of the ffow, and doesn't signiffcantly change the resulting 
length of the tail. This also translates into only a rela- 
tively small difference in the predicted observables, with 
changes of less than 30% in both X-ray and Ha lumi- 
nosities.^ The biggest difference is in the survival of H I 
clouds, with Tmin = 300 K predicting longer lived clouds. 
Still, the impact on H I observations is quite small. The 
ratios of the three observables is also quite constant, gen- 
erally to within 30%, as can be see from an inspection 
of the right two panels of Figures 1, 4, and 5, and from 
Tables 5 and 6. 

5.2. Resolution 

^ Note that this differs from the conclusion reached in TBIO, 
where we found that the Tmin = 300 K run predicted lower Ho; 
emission by more than an order of magnitude. This incorrect con- 
clusion was reached because of the cooling curve error for that run 
discussed in section 2. 



Resolution is most likely to affect the survival and 
structure of our dense clouds in the tail. The most direct 
results would be different amounts of neutral gas and Ho^ 
emission. A number of examinations of the survival of 
clouds with a variety of physics included have been per- 
formed. Mellema et al. (2002) and Yirak et al. (2009) 
include radiative cooling. Fragile et al. (2004) also dis- 
cuss how including self-gravity increases the lifetime of 
clouds. Nakamura et al. (2006) discuss the impact of us- 
ing smooth cloud boundaries on the growth of instabili- 
ties. They (and Yirak et al. 2009) find that a low density 
gradient results in slower growth of instabilities, which 
can retard cloud destruction. Our most dense clouds 
have an analytically calculated destruction timescale due 
to turbulent viscous stripping (using eq. (22) in Nulsen 
1982) of more than 1 Gyr, so resolution should not have 
a large impact over the timescales we consider in these 
simulations. We use the same cloud parameters as in the 
evaporation time calculation, although even with an or- 
der of magnitude lower density clouds (~ 10~^^ g cm~^), 
the destruction timescale is still about 200 Myr, which is 
much longer than the time at which we make our projec- 
tions. We use a relative velocity difference of 400 km s"-*^ 
(although we do not show velocity plots in this paper, 
see TBIO for a detailed discussion of the tail velocity). 
However, lower density clouds in our simulated tails cool 
and are compressed by the ICM, which may instead be 
destroyed by the ICM wind if we had better resolution 
that resulted in a steeper density gradient at the cloud 
edge. 

Although we do not perform a detailed resolution study 
of the runs in this paper, we do compare T3vl to a run 
with only 5 allowed levels of reffnement, to a minimum 
cell size of 76 pc, a factor of two worse in resolution. We 
compare outputs with the same amount of gas stripped 
from the galaxy. The total X-ray luminosity of the lower 
resolution tail is ~90% of the luminosity of T3vl (i.e. 
only a 10% change), while the total Ha ffux is only 28% 
of T3vl. The total HI ffux is also decreased by a sim- 
ilar amount (to 33% of the value in T3vl). For both 
resolutions, there is a good correspondence between the 
relative HI and Ha level and morphology; therefore, we 
argue that the ratio of HI to Ha is more robust than the 
absolute level of either. 

Both the small discrepancy in the X-ray luminosity 
and larger discrepancy in the Ha ffux can be explained 
by considering what gas is producing the emission. X- 
ray emission is produced by gas that is mixed with the 
ICM and is not localized near the dense clouds in the 
tail, and so is not very dependent on high resolution. 
One explanation for the slightly lower X-ray luminosity 
is that mixing happens more quickly and the gas is heated 
out of the X-ray bright regime. 

Unlike the diffuse nature of the X-ray emitting gas. 
Ha emission occurs at the edges of clouds. We ffnd that 
in projection the range of Ha intensities is the same as 
T3vl, meaning that the different resolution element size 
does not strongly affect the density and temperature at 
the edges of the dense clouds (and therefore does not 
strongly affect the Ha emission). The main difference 
is that there are fewer dense clouds in the less-resolved 
wind. In fact, the total H I column density in the low res- 
olution tail is about 30% of the total H I column density 
in T3vl-very similar to the Ha flux fraction. In projec- 
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tion, there are about half as many pixels with H I column 
densities over 10^^ cm~^. Since we expect the total H I 
column density to also depend on the number of clouds 
along the line of sight of the projection, we find that 
there should be roughly 26% as many clouds (by sim- 
ply squaring the 51% we find in the projection plane). 
Of course, clouds are more than a single pixel, so this 
is a rough approximation. However, this indicates that 
the amount of Ha emission closely follows the amount of 
H I column density, and likely the number of H I clouds. 
There are fewer H I clouds in the low resolution run be- 
cause with lower resolution the maximum density is lower 
and therefore the clouds are destroyed and mixed more 
quickly into the ICM. 

Deep, high resolution observations in H I will allow us 
to determine how many and for how long dense clouds 
survive, which will point to one of these mixing scenarios. 

6. CONCLUSIONS 

We have run detailed galaxy simulations including ra- 
diative cooling and made comparisons to the observed 
tail of ESO 137-001 to understand the physical mecha- 
nisms at work in the ICM. We compare three cases in 
which we vary our cooling floor between 8,000 K or 300 
K and the ICM parameters as shown in Table 3. Our 
main conclusions are as follows: 

1) The X-ray luminosity of our simulated tails are an 
excellent match to the X-ray luminosity of the tail of ESO 
137-001. This suggests that we are correctly modeling 
the phase distribution of gas in the tail, and that the 
mixms of the hotter stripped ISM (T > 10^ K) is beines: 
accurately modeled in the simulations. 

2) We find that bright X-ray emission depends upon 
a high surrounding ICM pressure, and find a minimum 
necessary pressure of 9 x 10~^^ erg cm~^. This conclu- 
sion agrees well with the local environment in clusters 
where bright X-ray tails are observed. 

3) We compare our Ha fluxes to the total diffuse Ha 
flux measure by Sun et al. (2007). As in our X-ray emis- 
sion, we find an excellent match between our simulations 
and observations. 

4) We predict that deeper observations will find H I 
gas to at least 40 kpc above the disk, because diffuse 
Ha emission directly traces neutral clouds (and has been 
observed to 40 kpc above the disk). The observations 
of molecular hydrogen by Sivanandam et al. (2009) 
strengthen this prediction. We also conclude that the 
mismatch between T3vl and Vollmer et al. (2001) in- 
dicates that the higher velocity cases better match the 
observations. 

5) Using the fact that Ha traces neutral gas and cal- 
culating the evaporation times of our simulated clouds 
(based on their sizes and densities), we calculated a max- 
imum efficiency for heat conduction of 24% (T3vh and 
T3vl) or 12% (T80vh). 

By modeling a simulation using a high ICM pressure, 
we have shown that X-ray emission can coexist with H I 
and Ha emission. We have also found that H I and Ha 
emission spatially coincide because Ha is mostly pro- 
duced at the edges of neutral clouds (Figure 3), while 
X-ray emission is generated in hot, diffuse gas that is 
mixing with the ICM. This is seen by the more even dis- 
tribution in the tail (Figures 1 and 3), and by the fact 
that the X-ray tail can be longer than the H I and Ha 



tails. 

Our excellent agreement with X-ray observations gives 
us confidence that we are correctly modeling the mixing, 
cooling and heating rate of X-ray emitting gas. This in 
turn means that heat conduction is not acting strongly to 
heat the diffuse gas in the tail, and that small scale turbu- 
lence (below our resolution scale) is not quickly mixing 
stripped gas. Both of these mechanisms would act to 
heat the gas in the tail out of the X-ray bright range and 
therefore lower the total luminosity, which would make 
our agreement to the observations worse. We cannot rule 
out the possibility that we have less gas mass at higher 
emissivities than in ESO 137-001. 

We also find excellent agreement between our simu- 
lated Ha flux and observations of diffuse Ha emission. 
We robustly conclude that diffuse Ha emission coincides 
with H I gas, and we find that Ha emission outlines H I 
clouds in all of our simulations. The agreement between 
the total Ha emission in our cases with different cooling 
floors underscores the fact that the edges of clouds, which 
are interacting with the ICM, emit the most strongly in 
Ha, not the central regions of the clouds, which are more 
likely to have radiatively cooled to the minimum allowed 
temperature. Although our agreement with observations 
indicates that we may also be correctly modeling the 
edges of dense clouds, we cannot dismiss the possibility 
that we are not fully resolving the cloud edges. In Section 
5.2 we discuss how lower resolution results in fewer H I 
clouds, and less Ha emission. Therefore, we cannot rule 
out that properties of our simulations that have not con- 
verged at our current resolution-the number of clouds, 
the rate of the decline of density at the edge of the cloud 
and the smallest scale of turbulent heating-have com- 
bined with the lack of heat conduction in our simulations 
to result in an accidental, and incorrect, agreement with 
observations. Observations of H I in the tail will provide 
an important check to these results. 

This work highlights the importance of comparing sim- 
ulations to detailed, mult i- wavelength observations of in- 
dividual systems. We are able to make predictions about 
this particular galaxy, such as the existence of H I gas in 
the tail and that its three-dimensional velocity relative 
to the ICM is probably larger than 1413 km s~^. How- 
ever, our simulation is not able to make any prediction 
about why there is a separated tail in ESO 137-001. We 
also do not reproduce the exact surface brightness distri- 
bution along the X-ray tail. As we have discussed, our 
goal was not to reproduce this specific tail. In order to 
do this, we would recommend modeling a smaller galaxy 
and matching the inclination angle between the galaxy 
and the ICM wind. 

Using our comparison with the observations of ESO 
137-001, we also draw more general conclusions about the 
importance of turbulence and the efficiency of heat con- 
duction in the ICM. We conclude that the mixing rate of 
the hot stripped ISM (T > 10^ K) is well-modeled in our 
simulations using only adiabatic compression and tur- 
bulent mixing down to the resolution of our simulations 
(38 pc). We call upon observers to test our predictions 
of where in clusters tails will be X-ray bright, and to use 
deep observations to verify the connection between Ha 
and H I gas. Observations of the H I tail of ESO 137-001 
can be used to test mixing of cold clouds into the ICM 
in our simulations. 
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